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Detecting the QCD phase transition in the next Galactic supernova neutrino burst 
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Predictions of the thermodynamic conditions for phase transitions at high baryon densities and 
large chemical potentials are currently uncertain and largely phenomenological. Neutrino observa- 
tions of core-collapse supernovae can be used to constrain the situation. Recent simulations of stellar 
core-collapse that include a description of quark matter, predict a sharp burst of Ve several hundred 
milliseconds after the prompt Ue neutronization burst. We study the observational signatures of 
that Ve burst at current neutrino detectors - IceCube and Super-Kamiokande. For a Galactic core- 
collapse supernova, we find that signatures of the QCD phase transition can be detected, regardless 
of the neutrino oscillation scenario. The detection would constitute strong evidence of a phase tran- 
sition in the stellar core, with implications for the equation of state at high matter density and the 
supernova explosion mechanism. 

PACS numbers: 14.60.Pq, 26.50.-f-x, 95.85.Ry, 97.60.Bw 



I. INTRODUCTION 



The majority of stars more massive than 8 Mq end 
their lives as core-collapse supernovae (SNe), with explo- 
sive kinetic energies ^ 10^^ erg. The explosion mecha- 
nism is most hkely related to the revival of the stalled 
bounce shock - the shock that forms when the collapsing 
iron core reaches nuclear densities, and stalls on its way 
out due to continuous energy loses via neutrino emission 
and dissociation of heavy nuclei. The detailed explosion 
mechanism has not been unambiguously identified yet. 
Given the surrounding matter envelope is opaque for pho- 
tons, neutrinos have been highly sought after for studying 
the physical conditions, dynamics of the collapse, and 
the SN mechanism. Indeed, current neutrino detectors 
are expecting high statistics from the next Galactic SN, 
and are on the verge of detecting the diffuse neutrino 
background from all past SNe [ll-y| • These would enable 
the explosion mechanism to be tested 043; '^^1^ 
neutrino mixing parameters to be constrained (see, e.g., 
Ref. 3). Moreover, neutrinos could potentially reveal 
new physics operating deep in the stellar core. 

Phase transitions have long been investigated in the 
context of SNe @ and neutron stars [13]. By includ- 
ing prescriptions of additional matter in the equation of 
state (EoS), many simulations predict the appearance of 
strange matter in the form of hyperons, a kaon conden- 
sate, or quark matter at supernuclear densities [l0l - [l2j . 
Bose condensates of pions have also been studied. In this 
paper we concentrate on recent developments in quark- 



hadron phase transitions. Initial attempts were based on 
general relativistic hydrodynamics using a parametrized 
equation of state @, discussing possible relations to the 
neutrino spectrum of SN1987A [l^. Applying a more 
sophisticated equation of state, the formation of a strong 
second shock wave was found as a direct consequence of 
the phase transition [l3]. However, these models could 
not predict the post-bounce neutrino signal due to the 
lack of neutrino transport. Phase transitions were also in- 
vestigated in the context of very massive (~ lOOM©) pro- 
genitor stars, where the time until black hole formation 
is shortened due to the softening of the equation of state 
during the quark- hadron phase transition [l2|, [H, ITg} . 
Recently, applying general relativistic radiation hydrody- 
namics based on three-flavor Boltzmann neutrino trans- 
port in spherically-symmetric simulations of low- and 
intermediate-mass progenitor stars, combined with the 
MIT bag model for the description of strange quark mat- 
ter, the formation of a strong second shock wave was 
confirmed [13 • The second shock accelerates at the sur- 
face of the protoneutron star and merges with the bounce 
shock, triggering an explosion, where otherwise no explo- 
sion could be obtained. As the second shock crosses the 
neutrinospheres, a second burst of neutrinos is released. 
Since the degeneracy is lifted in the newly shocked ma- 
terial and the electron fraction must be increased, the 
second burst is dominated by v^. 

Current neutrino detectors like IceCube and Super- 
Kamiokande are mostly sensitive to SN Pg's through in- 
verse beta decay reactions (p^p — > ne^) in the detector 
medium, while z/g's are detected by the subleading elas- 
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FIG. 1: Neutrino luminosities and mis neutrino energies as 
functions of time after bounce, sampled at 500 km radius in 
the comoving frame, for a 10 Mq progenitor star as modeled 
in [l3|: in solid (blue), i/s in dashed (red), and i^^/,- in 
dot-dashed (green). In contrast to the deleptonization burst 
just after bounce (f ~ 5 ms) the second burst at t ~ 257 — 261 
ms is associated with the QCD phase transition. The inset 
shows the second burst blown up. 



tic scattering channel. Encouragingly, the second burst 
shows up most prominently in Pg's, and its detection ap- 
pears feasible already. This contrasts with the delep- 
tonization burst which has been widely explored and pro- 
posed as a crucial tool to determine the flavor oscillation 
effects on the SN neutrino signal [H, [l^ ; since the delep- 
tonization burst consists of Ve^s, its detection requires 
larger neutrino experiments. If the Ve burst associated 
with a phase transition is observed, it would constitute a 
strong signature for the presence of a phase transition in 
the SN core, with important implications for new physics 
and the SNe explosion mechanism. We consider this 
intriguing possibility and investigate the capabilities of 
the largest existing neutrino detectors to unambiguously 
measure the Pe burst from a Galactic SN (assumed at a 
distance of 10 kpc unless stated otherwise). 



II. NUMERICAL MODEL 

We employ the neutrino emissions of Ref . [l3| , which 
simulated the core-collapse of massive stars using a quark 
matter equation of state based on the widely used MIT 
bag model. The main physical uncertainty in modeling 
the QCD phase transition is the critical density at which 
the transition sets in. In the simulations of Ref. [l^ , this 
is determined by the bag constant and the strange quark 
mass, chosen such that a quark-hadron transition sets 
in during the early post-bounce phase, close to nuclear 
saturation density. 

The neutrino emission of the reference simulation of 
Ref. [l3 is shown in Fig. [TJ The ordinary deleptonization 
burst is seen at i ~ 5 ms after bounce. At about the time 



when the quark matter core is created (t ~ 260 ms for 
the simulation shown) a strong second shock wave forms, 
which crosses the neutrinospheres in a few milliseconds 
and releases a second burst comprising of z^e's. The rise- 
time is a few ms, being related to the shock crossing 
the neutrinospheres. The duration is < 4 ms and the 
time-integrated energetics is approximately 5 x 10^^ erg 
in i^e's, ~ 1% of the total energetics. Additionally, the 
burst is accompanied by a sharp rise in the neutrino aver- 
age energies. The fluctuations in the neutrino luminosity 
shortly after the second burst are due to the appearance 
of an accretion shock when cooling at the neutrinospheres 
leads to the fallback of the innermost ejecta. The accre- 
tion rate and position of the accretion shock varies with 
time and modulates the luminosity until it settles to a 
quasi-stationary state on a time scale of about 100 ms. 

A second set of simulations was performed using a 
larger bag constant which results in a higher critical den- 
sity. The post-bounce time for the second burst is 448 
ms (see Table I of Ref. [13), but otherwise, the burst 
duration (^ 4 ms) and energetics 10^° erg in i/e's) 
are similar to those shown in Fig. [TJ The timing of the 
second burst after the deleptonization burst contains cor- 
related information about the quark-hadron equation of 
state, the critical conditions for the quark-hadron phase 
transition, and the progenitor model. 



III. NEUTRINO FLAVOR CONVERSIONS 

In order to determine the SN i?e signal observed at 
Earth, one must take into account flavor conversions oc- 
curring during propagation. In general, neutrino oscil- 
lation effects vary during the post-bounce evolution. In 
particular, collective flavor conversions in the SN [23| are 
forbidden during the second burst. From Fig. [TJ one re- 
alizes that during this phase the Ve flux is strongly en- 
hanced with respect to the non-electronic species and 
Dx (where x = M> ''') while the Ve flux is strongly sup- 
pressed. In this condition, the conservation of the lep- 
ton number prevents the collective I'eVe i^xi^x pair 
conversions [2l| . Therefore, only Mikheyev-Smirnov- 
Wolfenstein (MSW) flavor conversions occur while the 
neutrinos propagate through the stellar envelope [2^ . 

In inverted mass hierarchy (IH: < mi < TO2, where 
rrii is the neutrino mass), MSW matter effects in the SN 
envelope arc characterized in terms of the level-crossing 
probability Ph of antineutrinos, which is in general a 
function of the neutrino energy and of the 1-3 leptonic 
mixing angle ^13 [2^ . In the following, we consider two 
limits, namely Ph — when sin^ 6'i3 > 10^"^ (large) and 
Ph ~ 1 when sin^ 6*13 < 10~^ (small). Neglecting for 
simplicity Earth matter crossing effects, the electron an- 
tineutrino flux Fp^ at the Earth surface for IH with small 
013, is given in terms of the primary fluxes by (2^ 

Fp^ ~ cos2 ei2Fl + sin^ Oi^Fl , (1) 
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FIG. 2: Photon count rates at IceCube for the neutrino emis- 
sion shown in Fig. [T] The inset shows the second burst blown 
up, in tlie same axis units. 

where cos^ 9i2 — 2/3, with 6*12 being the "solar" mixing 
angle For IH with large 6*13 one has 

F,. Fl . (2) 

For normal mass hierarchy (NH: mi < 7712 < m^) the 
expression in Eq. ([T]) is applicable. The signal is expected 
to be larger for NH or IH with small 9i3 than that for IH 
with large 6*13, because F?^ > F^^. 

Outside the De burst, one finds F^^ > F^^ > F^^. This 
flux ordering allows for both collective and matter flavor 
transitions in IH. F[,^ is given by Eq. ([ij for large 613, 
and Eq. ^ for small ^13, i.e., effectively the case of large 
and small ^13 are exchanged relative to the case with 
no collective effects. For NH collective oscillations are 
usually unimportant, and we find the same result as in 
Eq. (P) for the final flux. 

IV. DETECTION OF THE QCD NEUTRINO 
BURST 

We analyze the detection of the QCD phase transi- 
tion induced i>e burst in two neutrino detectors, the km^ 
ice Cerenkov detector IceCube and the 32 kton water 
Cerenkov detector Super-Kamiokande. The detectors 
have distinct advantages - IceCube is greater in volume, 
while Super-Kamiokande allows spectral energy recon- 
struction. 



A. Detection at IceCube 

To estimate the event rate produced by inverse-/? decay 
in the ice, we follow the recent calculation of Halzcn and 
Raffelt Q. The complete detector will have 4800 optical 
modules and the data are read out in 1.6384 ms bins, 
implying a total event rate of 

R^^ = I860 bin-^Lss d^o (3) 



FIG. 3: Schematic representation of the relative time-ofTset 
between IceCube strings. 

where L53 is the Vf, luminosity in units of 10^"^ erg s~^, d 
is the distance to the SN in units of 10 kpc, and Ei^ is the 
i/g energy in units of 15 MeV. The spectrum-dependent 
factor in the expression is « 1.8 for the simulation in 
Ref . [17} . The signal rate has to be compared with the 
stochastic background rate Rq = 2.20 x lO'^ bin"-'^, with 
a rms fluctuation of 47 bin~^. The solid (blue) line in 
fig. [2] shows the event rate due to the De burst in IceCube 
for normal mass hierarchy [Eq. ([T])] in 2 ms bins. The 
dashed (green) line shows the expected background. In 
the narrow time window (^ 4 ms) where the burst would 
show up, one would collect 2 x 10* events, as opposed 
to ^ 8000 without the QCD phase transition. 

The vertical (red) band corresponds to the range of al- 
lowed signals for other possible oscillation scenarios. If we 
take the most pessimistic possibility that the oscillation 
effects are somehow arranged to produce the maximum 
possible signal outside the burst region and the minimum 
possible inside it, we get ^ 10* additional events in the 
burst bins over a steady background of ~ 10* SN events 
producing an excess with a significance of ~ lOOcr. The 
significance goes down at most quadratically with the SN 
distance, and even at 50 kpc (the LMC) remains signifi- 
cant. The signature of the Dg burst, associated with the 
QCD phase transition in the SN core, would be spectac- 
ular in IceCube. 

Since the temporal position of the j>e burst depends 
on the equation of state and the critical conditions for 
the phase transition, reconstructing the post-bounce time 
of the Ve burst is crucial. IceCube can reconstruct the 
bounce time to within ±3.5 ms at 95% CL for a SN at 
10 kpc 0, HI], hence the post-bounce time of the QCD 
burst may be measured to within about ±4 ms, which 
is sufficiently accurate to discriminate between different 
model predictions (see Table I of Ref. [T7|). 



Sub-Sampling at IceCube 

IceCube is usually quoted to have a time resolution of 
1.6384 ms. This is indeed the time-window over which 
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FIG. 4: Reconstructed neutrino event rates (note the subtle 
difference from Fig. [2]) at Super-Kamiokande, for the neutrino 
emission shown in Fig. [T] The inset shows the second burst 
blown up, in the same axis units. 



each optical module collects its data and transmits it to 
the data acquisition. The time synchronization of each 
module can however be controlled within ^ 3 ns relative 
to each other [25| . Thus the intrinsic timing of IceCube 
is much better. The coarser time resolution arises only 
due to the 1.6384 ms binning. This can be in principle 
remedied, if one staggers the signal collection windows. 
In Fig. |3l we schematically show the strategy. One could 
synchronize all modules on each string, but allow differ- 
ent strings to be offset in time by some fraction of the 
bin-size. For a sharp transient signal, different strings 
then attach a different time-stamp to the signal, and one 
would get a better time-resolution of the signal from a 
systematic comparison [26j . 

This procedure would not be very useful for weak sig- 
nals where one loses significance rapidly by sub-dividing 
the detector. However, the QCD burst signal is expected 
to be quite large, and this strategy will allow a precise 
measurement of the rise-time of the signal. 



B. Detection at Super-Kamiokande 

In the following paragraphs we investigate the detec- 
tion at Super-Kamiokande. The simplest procedure to 
perform is to similarly time-bin the events. We use bins 
of 5 ms, which is comparable to the duration of the sec- 
ond burst. Fig. |4] shows the reconstructed neutrino event 
rates expected in Super-Kamiokande for the normal mass 
hierarchy in square points (blue) with statistical error 
bars. The range of allowed rates for other possible os- 
cillation scenarios is shown by the shaded band (red). 
The QCD burst would produce ~ 60 excess events in the 
two relevant bins. The statistics are encouraging enough 
that the QCD-burst would be identifiable even at a SN 
distance of 20 kpc, which would contain 90% of the 
Galactic SNe [23|. However, the smaller the signal to 
noise ratio, the larger the risk that fluctuations in the ac- 
cretion flow, which also produce variations in the electron 



flavor neutrino luminosities (see, e.g., [28|) can obscure a 
clear signal of the phase transition. 

Super-Kamiokande's energy resolution can reveal a 
second distinguishing feature of the QCD phase transi- 
tion burst - the increase of the average neutrino ener- 
gies. In particular, the largest energy increase is seen in 
the Dx, which rapidly increases from 20 MeV to 35 MeV 
(sec Fig. [1} . For the inverted mass hierarchy with large 
013, this energy increase would appear in the observed 
De- Super-Kamiokande's energy resolution for positrons 
in this energy range is 10% or better The De en- 

ergy can be determined from the positron energy, and 
hence the energy resolution of Super-Kamiokande per- 
mits the detection of this secondary effect due to the 
sudden compactification of the progenitor core during 
the quark-hadron phase transition. A third feature of 
the QCD phase transition is the transient increase in 
the Pe energies and luminosity, occurring ~ 0.1 s after 
the phase transition, as a result of accretion onto the 
remnant. This feature is most readily observable for the 
normal hierarchy. 



V. DISCUSSION AND SUMMARY 

Core-collapse SNe can be used to probe the state of 
matter at high baryon densities and finite temperatures, 
beyond the reach of heavy-ion collision experiments at 
RHIC and the LHC, and complementary to those that 
will be studied with FAIR at the GSI Darmstadt. Using 
the MIT bag model, Ref. [l3 investigated the quark- 
hadron phase transition during the early post-bounce 
phase of core-collapse SNe, demonstrating that it would 
cause a prominent burst of De- We showed that this burst 
from the next Galactic core-collapse would be clearly de- 
tectable at IceCube and Super-Kamiokande, regardless 
of neutrino mixing scenarios. 

We have discussed three observational features of the 
QCD related i>e burst. The phase transition results in 
a clear and sharp rise in the signal rates at neutrino 
detectors, with Icc-Cubc providing accurate measure- 
ments of the luminosity and timing of the burst. Super- 
Kamiokande also measures the energies of the neutrinos, 
revealing the increase in average neutrino energies pre- 
dicted to occur during the burst. Finally, fluctuations in 
the accretion rate after the phase transition are also de- 
tectable. The accurately measured energetics and timing 
would reveal information about the critical conditions for 
phase transitions in the protoneutron star. 

However, other transitions may occur in the extreme 
densities of the neutron star, e.g. pion and kaon conden- 
sations [30j . Indeed, any phase transition taking place in 
the early protoneutron star phase, through an extended 
mixed phase, rnay produce signals similar to those stud- 
ied by Ref. [I4I and this paper. Therefore, the detection 
of a burst of during the early phase of a SN would 
provide a strong signature of a phase transition, and con- 
strain the required conditions for its onset. The energet- 
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ics and timing of the i/g burst would help obscrvationally 
constrain the conditions required for the phase transi- 
tion. In this regard we note that the physics of the phase 
transitions and whether a sharp burst can arise therein 
remain to be studied systematically. 

In addition to SN physics, the short rise time of the 
burst, combined with high statistics and energy measure- 
ment capabilities of SupcrKamiokande, would allow the 
neutrino mass to be limited by the energy-dependent de- 
lays in arrival times. A massive neutrino with mass m 
and energy traveling a distance D, will experience a 
delay relative to a massless neutrino, of 



At K. 0.515 



(m/eV)2(i:)/10kpc) 
(£;/MeV)2 



■ s. 



(4) 



Considering the detector response and statistical errors, 
the mass limit from the QCD burst would be < 1 eV, 
assuming a rise time of 1 ms and that the rise includes 
100 events [3l|. This limit would be comparable with 
the most optimistic forecast obtained through a future 
measurement of SN neutrinos in SupcrKamiokande [s^] ■ 
However, the detection of the first few events plays a key 
role in obtaining a strong limit, implying a large fluctua- 
tion on the mass bound. The mass limit achievable with 
the QCD burst would also be competitive with the less 
direct cosmological limits on neutrino masses [33| . 

This antineutrino burst signal also has implications for 
coincidence measurements with gravity wave detectors. 
The exceptional timing provided by the narrow QCD 



burst will improve comparison with gravit y w ave detec- 
tors, for the procedure suggested by [H,!!!]- A clear 
separation of the different bounce, accretion, and cooling 
phases of the explosion signal will allow for more sophis- 
ticated analysis of the signal. 

In conclusion, our study confirms the crucial role of as- 
trophysical multi-messengers played by neutrinos during 
a stellar collapse. In particular, it is fascinating to real- 
ize that weakly-interacting particles, like neutrinos, could 
become a powerful tool to probe the nature of strong in- 
teractions in extreme conditions. 
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